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(57) Abstract 

A non-invasive system and procedure for 
deriving the blood gas content fox a patient. The 
system measures the carbon dioxide concentra- 
tion of the expiratory breath relative to volume. 
This data is then processed to* derive arterial 
blood gas levels of carbon dioxide. If data sam- 
pling is in the time domain, the processing shifts 
the data from the time domain to the volume do- 
main. The processing also iteratively assesses 
the significance of numerous variables. The re- 
sulting relationship provides a fast and accurate 
measure of blood gas content for both healthy 
and diseased lung patients. 
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WO 96/24285 

NON-INVASIVE ESTIMATION OF ARTERIAL BLOOD* GASES 

BACKGROUND OF THF- INVENTION 
5 Field of the Invention : The present invention relates to methods and 

apparatus for assessing the ventilatory status of a patient. More particularly, the 
present invention provides a system for implementing a non-invasive procedure for 
estimating the amount or concentration of dissolved carbon dioxide within the 
arterial portion of the vasculature. The arterial carbon dioxide content, expressed 

10 as a partial pressure, i.e., pC0 2 is an important measure of ventilatory status which 
ultimately reflects pulmonary health. 

State of the Art : Physicians and other health care providers often use 
elevated arterial pC0 2 (PaC0 2 ) as an indicator of incipient respiratory failure. In 
this regard, the determination of PaCO,, is useful in optimizing the settings on 

15 ventilators and detecting life-threatening blood gas changes in an anesthetized patient 
undergoing surgery. The traditional method for obtaining arterial blood gas values 
is to extract a sample of arterial blood and measure the partial pressure of carbon 
dioxide using a blood gas analyzer (PaC0 2 ABC ) . Arterial puncture has inherent 
limitations: 1) arterial puncture carries a degree of patient discomfort and risk, 2) 

20 handling of the blood is a potential health hazard to health care providers, 3) 
significant delays are often encountered before results are obtained and, 4) 
measurements can only be made .intermittently. 

Continuous invasive monitoring requires in-dwelling arterial lines which 
entail inherent problems. These include sepsis, slow response times, and signal 

25 decay. The nature of this monitoring systems excludes its use under routine care 
and is generally restricted to intensive care units within a hospital facility. 

There have been attempts to assess PaCO = levels indirectly, including a 
technique known as capnography. The approach utilized in capnography involves 
tracking patient exhalation and measuring expiratory gas CO; concentration against 

30 time during one or more respiratory cycles. The resulting relationship is plotted to 
create a graph depicting three distinct phases in breath C0 2 gas concentration during 
the patient exhale cycle. (See, Figure 1). Typically, the three phases reflect the 
clearing of the conducting airways which do not normally participate in gas 
exchange (i.e., airway dead space) (phase I) followed by the exhalation of air from 
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conducting airways dynamically mixed with lung gases from the active (alveoli) 
membrane surfaces within the lung that have undergone gas exchange with arterial 
blood (Phase II). The final phase (phase HI) reflects the exhalation of unmixed gas 
from regions of the lung which are normally in active exchange with the alveoli 
5 tissue. Phase m thus closely resembles (in healthy patients) gas properties 

associated with arterial blood in contact with the lung for gas exchange (CO, release 
and 0 2 absoiption). In normal lungs, the graph line of Phase III is substantially 
level (slope = 0) since ventilated and perfused alveolar regions are closely matched. 
In a diseased lung, the Phase m line may not appear level due to a mismatch in 
10 ventilation and perfusion of this lung region. Table I below: 

TABLE I 

Phase I Represents C0 2 -free gas expired from the airway conduction 

structures where gas exchange does not occur, 

Phase D The S -shaped upswing in C0 2 concentration (expressed as a 
15 percent) represents the transition from airway to alveolar gas, 

and 

Phase IE The alveolar plateau representing C0 2 rich gas from the 
alveoli. 

In the past, capnography has utilized the peak or end-tidal (PetCOa) values as 
20 an estimate of PaC0 2 . PetC0 2 is a measure of the mean alveolar partial pressure of 
carbon dioxide from all functional gas exchange units of the lung. PetCQ. obtained 
from capnography is a measure of mean alveolar pC0 2 , which value approximates 
PaC0 2 in normal lungs. Because C0 2 readily diffuses across the alveolar-capillary 
membrane, the PetC0 2 closely approximates the PaC0 2 with normal ventilation- 
25 perfusion. The difference between PetC0 2 and PaCo, is primarily a function of the 
proportion of the lung where gas exchange does not occur (Fletcher, R., Johnson, 
G., and Brew, J., The Concept of Deadspace with Special Reference to Single 
Breath Test for Carbon Dioxide. Br. J. Anaesth. . 53, 77, 1981). In patients 
afflicted with a lung disease, there often exists a proportional increase in the region 
30 of the lungs where gas exchange does not occur. This increase in so-called 

"alveolar dead space" results in a significant difference between peak C0 2 (PetC0 2 ) 
obtained from capnography and elevated arterial C0 2 (PaC0 2 ). 
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Other techniques have been utilized for assessing patient blood gas levels 
with mixed results. Transcutaneous sensors measure tissue pCOj diffused through 
the heated skin but have practical and theoretical limitations. Oximetry is a widely 
used, non-invasive method for estimating the arterial oxygen carried on hemoglobin. 
5 For example, U.S. Patent Nos. 4,759,369, 4,869,254 and 5,190,038 describe pulse 
oximeters which measure the percentage of hemoglobin which is oxygenated. 
However, neither of the aforementioned techniques measures the amount of 
dissolved oxygen present, nor the amount of oxygen carried when hemoglobin levels 
are reduced. Low hemoglobin levels are found when there is a significant blood 

10 loss or when there is insufficient red blood cell formation. Additionally, oximeter 
readings are specific to the point of attachment, which is typically the finger tip or 
ear lobe, and may not reflect the oxygen level of vital organs during conditions such 
as shock or hypothermia. 

There remains a significant need in the art for an accurate, non-invasive, 

15 sensitive method for accurately determining the levels of arterial blood gases. As 
will be seen hereinafter, the instant invention sets forth a non-invasive system to 
overcome the problems of the prior art. 

SUMMARY 

20 The present invention provides a system to rapidly and accurately derive a 

patient's arterial carbon dioxide concentration and employs a non-invasive method 
for monitoring arterial partial pressure of carbon dioxide in a patient as an indicator 
of ventilatory status. 

The present invention also comprises a system for detecting expiratory C0 2 

25 concentration and volumetric rate data and accurately deriving actual arterial pC0 2 
based thereon. The system affords non-invasive, substantially real time, 
determination of blood gas concentrations as derived from current expiratory data as 
correlated with processed data collected from past expiratory measurements. 

The present invention converts expiratory data from a time domain to a 

30 volume domain or by accumulating data in the volume domain in the first instance. 
The arterial C0 2 partial pressures of a patient can then be ascertained by selectively 
analyzing the slope and intercept values associated with Phase II and Phase in 
expiratory data in the volume domain. If time domain data is converted to the 
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volume domain, the conversion may be said to "normalize" the data by placing it in 
the volume domain, wherein time-dependent respiratory differences between patients 
are eliminated and a standardized gas concentration-to-incremental breath volume 
relationship is achieved. 
5 The present invention also provides for the accurate determination of arterial 

C0 2 partial pressures by measuring expiratory gas data and statistically filtering this 
data to ascertain readings having the highest correlation to actual pulmonary 
performance with a statistically significant level of confidence. 

The above and other advantages of the present invention are realized in a 

10 specifically delineated gas analysis and data processing system operated in 

accordance with select data qualifying and enhancing procedures. In particular, the 
inventive system provides for the collection of concise expiratory data from a 
patient undergoing treatment. This data includes details on C0 2 gas partial 
pressure, concentration, and total gas volume sampled during breath exhaust cycle. 

15 Multiple readings are made to enhance accuracy. The expiratory data is converting 
to a volume domain from the time domain. The expiratory data is then charted to 
establish the three aforementioned distinct phases within the expiratory cycle. The 
associated linear details of these three phases are extracted and used to project 
current arterial C0 2 partial pressure, PaC0 2 . This value is then utilized to quantify 

20 pulmonary performance and/or determine the existence of lung failure or distress. 

In accordance with the varying aspects of the present invention, the system 
includes a sophisticated artificial intelligence engine that iteratively analyzes many 
separate and distinct measurements and the calculated values associated with the 
expiratory data. Based on these permutations, the engine quantifies those measured 

25 and derived values having the highest correlation to the actual arterial C0 2 partial 
pressures, resulting in a fixed relationship including specifically weighted variables 
for projecting arterial C0 2 partial pressures. This relationship, expressed as a 
vector, is implemented according to system parameters in actual patient monitoring 
during surgical procedures and other periods of time associated with potential 

30 pulmonary failure. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
The foregoing features of the present invention are more fully and readily 
understood from the following detailed description of a specific illustrative 
embodiment thereof, presented hereinbelow in conjunction with the accompanying 
5 drawings of which: 

Figure 1 graphically illustrates three phases of the C0 2 /volume curve in a 
healthy person; 

Figure 2 represents the normal expiration of COj plotted against the expired 
volume and the represented dead spaces; 
10 Figure 3 is a schematic illustration of a system in accordance with a 

preferred embodiment of the present invention; 

Figure 4 is a flow chart of a preferred method of the present invention; 

Figure 5 is a capnograph for a healthy individual; and 

Figure 6 is a capnograph of a person suffering from lung distress. 

15 

DETAILED DESCRIPTION OF THE PREFE PPFD EM BODIMENT 
In overview, the present invention comprises a non- invasive system and 
procedure for deriving the gas content of arterial blood through analysis of the 
expiratory gas concentrations and volumes. The expiratory gas concentrations and 

20 volumes may be collected in a time domain, typically using a capnometer and 

pneumotachometer. Alternatively, the data may be collected directly in the volume 
domain with a time delay or phase-shift correction between the pneumotachometer 
and capnometer data. The raw data for both pC0 2 and volume are digitized at a 
frequency high enough to avoid aliasing. In adults with normal ventilatory 

25 frequencies, the data is typically collected at 100 Hz, producing a data point every 
0.01 seconds. If taken in the time domain, the measurements are converted from 
the time domain (typically 1 point per 0.01 seconds) to a volume domain (typically 
1 point per 1 ml of expired volume) utilizing a polynomial fit. 

In the volume domain, the slope of the Phase II curve segment is calculated. 

30 The initial point of Phase II is where the curve transcends a threshold value 
(typically 0.5% for adults). The final point of Phase II is defined, in this 
application, as the point at which the curve deviates from linearity by a specified 
amount (typically 5% for adults). The threshold and deviation values are based on 
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lung size and respiration rate. The slope of the Phase m segment is calculated, in 
an analogous manner, by calculating from the last data point of expiration toward 
Phase II until the curve deviates from linearity by a fixed amount (typically 5% for 
adults). The slopes of the Phase II and Phase ID curve segments are used as initial 
5 input variables, which variables are combined into a vector with other parameters to 
produce an estimate of content of the gas in the anerial blood, i.e., the partial 
pressure of C0 2 . or PaC0 2 . 

The following abbreviations and definitions are salient to invention 
description: 

10 ABBREVIATIONS AND DEFINITIONS 

C0 2 Carbon dioxide 

pC0 2 Partial pressure of carbon dioxide 

PetC0 2 End tidal C0 2 

PEC0 2 Expiratory C0 2 

15 PaC0 2 Arterial partial pressure of C0 2 

PaC0 2 ABG Arterial partial pressure of C0 2 measured by an arterial 

blood gas analyzer 

%C0 2 Carbon dioxide value expressed as a percentage of the 
total gas content 

20 COPD Chronic obstruction pulmonary disease 

ml Milliliters 

fds Fowler dead space (Airway dead space) 

12 Intercept of Phase II on C0 2 axis 

13 Intercept of Phase m on C0 2 axis 
25 m2 Slope of Phase II 

m3 Slope of Phase m 

ph2i Phase II deviation index 

ph3i Phase m deviation index 

Tb Phase transition bend 

30 Tbi Phase transition bend index 

TV Tidal volume of air exhaled 

ang Angle between Phases II and III 

As noted above, expired C0 2 pattern obtained from capnography may be 

graphically depicted as a three-phase curve. In a normal person (Figure 5), the CO : 

35 curve forms a plateau at the start of Phase m and reaches a value (PetC0 2 ) 

approximately equal to the PaCOj at the end of Phase m. In a patient with 

increased volume of the lung where gas exchange does not occur (physiologic or 

alveolar dead space, Figure 6), as typically occurs in COPD, the curve does not 

plateau and the final PetC0 2 values are a less reliable indication of the actual 



WO 96/24285 



PCT/US96/01538 



PaC0 2 . The instant invention determines an estimated value which closely 
approximates the actual PaC0 2 under these and other conditions. 

The inventive system is illustrated in Figure 3 and comprises a differential 
pressure flowmeter or pneumotachometer (item 15), a flow signal unit (item 14), a 
5 gas sensor (item 16), a C0 2 signal unit (item 13), a processor (item 12) and a data 
display (item 11). This system can be used with or without mechanical ventilation 
of the patient. 

Many devices for measuring the volume of a person's expiratory breath 
already possess the capability to integrate a measured flow and can be used in this 
10 invention- Typically, flow -measuring devices use one of the following methods to 
determine flow: 

1 . measurement of pressure drop or differential pressure across a fixed 
resistance (differential pressure flowmeter or pneumotachometer), 

2. measurement of the temperature change of a heated wire cooled by 
15 the airflow (hot wire anemometer), 

3. measurement of frequency shift of an ultrasonic beam passed through 
the airstream (ultrasonic Doppler), 

4. counting the number of vortices shed as air flown past a strut (vortex 
shedding), or 

20 5. measurement of transmission time of a sound or heat impulse created 

upstream to a downstream sensor (time of flight device). 

Alternatively, volume may be measured directly by counting revolutions of a 

vane placed in the flow path (spinning vane). A discussion of the aforementioned 

devices and associated technology can be found in Sullivan, et aL, 

25 Pneumotachographs: Theory and Clinical Application . Respiratory Care, Vol. 29-7, 
pages 736-749 (1984), which is incorporated by reference herein. Examples of 
known differential pressure flowmeters include those described in U.S. Patent Nos. 
4,047,521, 4,403,514, 5,038,773, 5,088,332, 5,347,843 and 5,379,650, the 
teachings of which are incorporated by reference herein. 

30 The exemplary device for respiratory flow measurement is the differential 

pressure flowmeter or "pneumotachometer" (Figure 3, item 15), which provides a 
pressure differential indicative of respiratory flow, the differential being converted 
via transducers in flow signal unit (item 14) to electrical signals representative of 
the relationship between respiratory flow and pressure differential. The flowmeter 

35 (item 15) is manufactured and sold by Novametrix Medical Systems, Inc., 



BNSDOCID: <WO 9624285A1_I_> 



WO 96/24285 PCTAJS96/01538 

Wallingford, Connecticut. However, any of the aforementioned types of flow 
measurement devices may be utilized in the inventive system. 

Sensors capable of measuring carbon dioxide content in a person's expiratory 
breath are well known: The currently preferred device for measuring carbon 
5 dioxide content is a gas analyzer of the type employing non-dispersive infrared 
radiation which presents data representing the %C0 2 (or pC0 2 ) of a sample of 
exhaled breath. Examples of known infrared gas analyzers include those described 
in U.S. Patent Nos. 4,859,858, 4,859,859, 4,914,720, 4,958,075, 5,146,092, 
5,153,436, 5,206,511 and 5,251,121, the teachings of which patents are 
10 incorporated by reference herein. Other technologies used to measure the 

concentration of carbon dioxide such as Raman spectroscopy and mass spectroscopy 
can also be used in the present invention. 

The exemplary gas sensor (Figure 3, item 16) capable of measuring carbon 
dioxide content in a patient's exhaled breath is available from Novametrix Medical 
15 Systems, Inc., Wallingford, Connecticut, under the trade name CAPNOSTAT™. 

Other methods of measuring carbon dioxide content both at the airway (mainstream) 
or by removing a sample (sidestream) may be used in the present invention. 

Gas analyzers as described above employ non-dispersive infrared radiation to 
measure the concentration of a selected gas in a mixture of gases. The infrared 
20 radiation can be emitted from a thick film source and focused by a mirror to pass 
through the mixture of gases being analyzed. After passing through the gas 
mixture, the infrared beam is passed through a filter which reflects all of the 
radiation except for those in the narrow bands centered on a wavelength which is 
absorbed by the gas of concern (such as COj) in the mixture being analyzed (such 
25 as the air flow from a person's expired breath). This narrow band of radiation, 

which typically extends approximately 190 angstroms to each side of the wavelength 
on which the radiation is centered, is allowed to reach a detector which is capable 
of producing aiT electrical output signal inversely proportional to the magnitude of 
the infrared radiation impinging upon it, as the radiation in that band is attenuated 
30 to an extent which is proportional to the concentration of the designated gas in the 
mixture of gases being analyzed. The strength of the signal generated by the 
detector is consequently inversely proportional to the concentration of the designated 
gas and can be inverted to provide a signal indicative of that concentration. The 
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processor can be either a personal computer with a suitable digital interface for 
receiving the digital signals from the C0 2 signal unit (Figure 3, item 13) and flow 
signal unit (Figure 3, item 14), or any specially designed processor capable of 
calculating the vectors as disclosed further herein. 
5 As discussed above, capnography systems in the prior art have attempted to 

estimate the PaC0 2 by using the peak or end-tidal value for a single breath 
(PetC0 2 ). Such systems are effective at tracking gas changes in normal lungs but 
have been shown to be unreliable in diseased lungs such an COPD or when a 
significant pulmonary shunt exists. The inadequacy of PetC0 2 to measured PaCO : 

10 is in part attributed to regions within the lung with high ventilation to perfusion 
ratios. Exhaled gases from regions where gas exchange does not occur due to 
inadequate perfusion reduce the obtained PetC0 2 . For example, in the diseased 
lung, the increased slope of Phase lH due to impaired gas exchange makes the 
PetC0 2 valve a less reliable indicator of PaC0 2 . 5 

15 In the instant invention, the digitized expired C0 2 and flow data, if taken in 

the time domain, is converted to the volume domain to account for variations 
between different people and thereby improve accuracy. Alternatively, flow signal 
unit 14 may be employed to integrate a flow rate signal on a continuing basis into 
volume, C0 2 readings then being taken at predetermined volumetric intervals, the 

20 phase shift between the two sensors (including signal processing time) being 

corrected for as data is taken or subsequently processed. Alternatively, all data may 
be taken in the time domain and converted to the volume domain after the fact by 
processing unit (Figure 3, item 12). Finally, as previously noted, flow volume may 
be measured directly, as by a spinning vane device. C0 2 concentrations versus 

25 respiratory flow or volume may be depicted as a curve including a series of units, 
each defined by a pair of points. In the time domain, the progression from unit to 
unit is based on fractions of a second, regardless of the quantity of air expelled 
during that time period. The effective rate of sampling depends on the patient *s rate 
of respiration. In the volume domain, however, each unit to unit movement is 

30 based on a unit of volume of air expelled, regardless of the expiratory time. By 
employing a volume domain or by transferring the data from a time domain to a 
volume domain a more physiologic relationship is expressed. Toward the end of 
the breath, relative percentage increments in expiratory time and volume vary 
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of the breath receives a weighting more indicative of its physiologic importance 
regardless of the health of the patient. For example, in a mechanically ventilated, 
chemically paralyzed, patient with COPD experiencing bronchospasm, the 
exhalation period is prolonged due to the narrowing of the airways (Figure 6). 
5 With pharmaceutical intervention, expiratory resistance can be reduced, thereby 
reducing time required to complete exhalation. The exhalation time greatly varies 
between the non-medicated and medicated patient, yet the exhaled C0 2 to total gas 
volume exhaled ratio is relatively constant. Plotted in the volume domain, the 
capnograph of non-medicated and medicated patients would be the same. 

10 Conversely, in the prior art systems (time domain), measurements taken before and 
after medication would vary greatly. 

As shown in Figure 3, as a person exhales, the C0 2 sensor (16) measures 
the pCQ 2 in the person's expiratory breath. Virtually simultaneously, the 
pneumotachometer (15) measures the flow of the person's expiratory breath using 

15 differential pressures across a fixed resistance as previously referenced. The 
differential pressure values from the pneumotachometer (15) are converted to 
electrical signals and digitized in flow signal unit (14) and the analog signals from 
C0 2 sensor (16) are digitized in C0 2 signal unit (13). The digital signals are 
processed, as further disclosed herein, in the processor (12) and displayed via 

20 printer, VDT, LED or other display devices as known in the art (11). The 

estimated PaC0 2 value is displayed after several breaths. The actual lag time is 
dependent on the consistency of the data with a minimum of six breaths initially and 
every third breath thereafter. This data flow is illustrated, as part of the complete 
conversion routine, in Figure 4. As stated heretofore, the COj sensor (16) and 

25 associated signal unit (13) measures the C0 2 content of the patient's breath. The 
output from the flow sensor (15) is representative of the flow of the expiratory 
breath. The pneumotachometer (15) outputs are converted via transducers in flow 
signal unit (14) to electrical signals which are digitized as time dependent signals 
representing the flow rate of the patient's breath, and subsequently integrated into 

30 volume. These signals are sent to the processor (12). 

As shown in the flow chart in Figure 4, the flow data and C0 2 data are 
isolated breath-by-breath and converted to the volume domain in order to express a 
more physiologic relationship. The volume domain progression from unit to unit is 
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based on a unit of volume expired without regard to time. As a preliminary 
programming step, a variety of variables are computed which depict features of 
Phase n, the transition period, and Phase m. The computed variables are selected 
to include all features and idiosyncrasies of the capnograph. The derived variables 
5 are of two types: 1) those with established physiologic importance, and 2) those 
which mathematically represent features of the capnograph whose physiologic 
significance is less clear. These derived variables have utility in normalization of 
the data to allow for diverse lung sizes and respiratory patterns. Artificial 
intelligence is then used to evaluate each variable and assign a mathematical weight. 
10 In particular, a neural network is employed to insure the variables are evaluated 
without imposing an initial bias. These variables are computed as described 
hereafter. 

Suitable software systems are readily available in the marketplace, and are 
exemplified by the neural network presently offered by NeuralWare, Inc. located at 

15 Penn Center West, Building IV, Suite 227, Pittsburgh, Pennsylvania 15276. 

Computation of Intermediate Variables 
A threshold level of C0 2 is detected by finding the point were the curve 
transcends from a value below 0.5% to a value above 0.5%. In the instant 
invention, this point is used as the initial point of Phase n. From the threshold 

20 point, subsequent C0 2 data points are adjoined and tested for linearity to the data 
point where the signal deviates from a linear path. The amount of allowed 
deviation is typically 5% in adults. The point where the deviation occurs marks the 
termination of Phase II and the slope of the Phase II segment is derived. As 
described further herein, the Phase II slope is used later. The slope of Phase m is 

25 determined in an analogous manner by starting at the last data point of expiration 
(Pet COj) and regressing toward the termination point of Phase II. Once the Phase 
n and ID slopes are computed, the remaining input variables are derived as the 
following: 

TARTF.n 

The Fowler dead space volume is determined by examining 
the volume that gives equal area between the start of the Phase 
II data and the start of the Phase HI data. This computation is 
done by fractional difference. (Fowler W.S. Lung Function 
Studies n. The Respiratory Dead Space. Am. J. Phvsiol. 
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154: 405, 1948 - the teachings of which are hereby 
incorporated by reference). 

The C0 2 axis intercept of Phase IL The effective Phase II 
slope (mO) is used at the Fds volume. The equation is: i2 = 
modeIC0 2 (Fds) - mO * Fds where modelCOj (Fds) is the 
model evaluated at Fds. 

The slope of Phase m. 

The C0 2 axis intercept of Phase HI. The Phase EI slope is 
used at the end tidal volume. The equation is: i2 = 
modelC0 2 (PetC0 2 ) - m3 * PetCC^ where modelC0 2 (PetC0 2 ) 
is the model evaluated at PetC0 2 . 

The angle in degrees between the Phase II and Phase IH. The 
effective slope of Phase II (mO) is used. The angle between 
the lines is: (tan(m3) - tan(mO) + x) * (180/r). 

The end tidal value of pC0 2 . The last element in the C0 2 (vol) 
data array is used for PetCo 2 . 

The total volume of air exhaled expressed in ml. 

The Phase II deviation index. This is smallest volume at 
which the Phase II line deviates from the smoothed C0 2 data. 
The comparison is done by fractional difference. 

The Phase IH deviation index. This is the smallest volume at 
which the Phase HI line deviates from the smoothed C0 2 data 
by the phase transition width (Tw). The negative of the 
second derivative of the model curve exhibits a peak in the 
region between Phase II and Phase m. The width of this peak 
is Tw and is determined as the full width at half the maximum 
of the peak. All second derivative curves have peaks that fall 
below half the maximum of increasing volume. However, the 
peaks in some derivative curves do not fall to 1/2 the 
maximum for decreasing volume. Thus, Tw is found from the 
right side half-width of the peak at half-maximum and 
multiplying by 2 to construct the full-width. 

The phase transition bend. This value is the maximum of the 
negative of the second derivative of the model curve. It 
reflects the measure of the sharpness of the transition. 

The phase transition bend index. This value is the index 
(volume) at which the phase transition bend (Tb) occurs. 
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ModAng A modified version of the angle in degrees between the Phase * 
n and Phase in lines. The modification scales the actual 
angle to enhance differences in the model fit curves. The 
scaling is heuristic. The modified angle is computed as: 
5 (tan(400 * m3) - tan(400 * mO) + x) * (180/ x). 

The product of these intermediate variables and their assigned mathematical 

weights are used in arguments in a hyperbolic scaling function to compute the final 

PaC0 2 . The mathematical weights for each variable were initially defined from a 

set of 100 derived variables used as inputs to a fully connected, back-propagation, 

10 neural network and the actual arterial carbon dioxide used as the desired output. 

Once trained for 75,000 iterations, the resultant weights were extracted and used for 

subsequent derivations where the actual PaC0 2 was to be determined. 

System operation based on the stored vector arguments are shown in the 

flow chart in Figure 4. As described herein, the system inputs the flow data and 

15 C0 2 data as isolated breath-by-breath and converted to the volume domain in order 

to express a more physiologic relationship. The volume domain progression from 

unit to unit is based on a unit of volume expired without regard to time. Referring 

to Figure 4, logic begins at start block 100, and begins a processing loop defined by 

block 110, defining a sequence of 10 iterations (exemplary). At block 120, the 

20 digitized respiratory data is loaded as variable RESP(I), converted to the volume 

domain at block 130, and then mapped, setting forth the three aforementioned 

phases of the capnograph, block 140, VRESP(I). 

At block 150, Phase n threshold is detected, THRESH(I), followed by block 

1 60 wherein the detection of the Phase n deviation to Phase m, DEVIA(I) is 

25 effected. Based thereon, the variables for the capnograph are determined, VAR(I), 

at block 170. The stored arguments based on the assessed weighting for the vector 

are recalled at block 180, and the final vector generated, VEG(I). Based thereon, 

the PaC0 2 is derived and stored for the Ith iteration, block 190. At block 200, this 

process is repeated for the next sequence of respiratory data from the patient, for 10 

30 iterations or breath cycles (exemplary). 

Starting at block 210, the data is refined by discarding four of the ten 

readings (two high and two low), XPaCO^, and the resulting six iteratives are 

statistically assessed and accepted if meeting a minimum confidential level (95 or 

99%), by way of example. If accepted, the mean value of the remaining iterations 
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is displayed, MPaCO a (I), block 220. The process continues giving pseudo real time 
data, as noted at block 230. 

The results of system operation is depicted for the following patients: 

Table JU 

5 PaCOj 



10 



15 



20 



25 



30 



Derived PaCOj 


PaC0 2 ABG 


mean 




err 


Patient No. 


mm Hg 


mm Hg 


std deviation 


mm Hg 


1. 


33.00 


34.30 


2.40 


1.30 


2. 


33.00 


31.33 


1.84 


1.67 


3. 


33.00 


32.20 


2.33 


0.80 


4. 


33.50 


32.00 


1.77 


1.50 


5. 


38.30 


39.67 


1.36 


1.39 


6. 


40.00 


38.37 


0.86 


1.63 


7. 


40.00. 


41.38 


0.62 


1.38 


8. 


42.00 


42.63 


0.15 


0.63 


9. 


42.00 


44.01 


1.60 


2.01 


10. 


44.00 


45.47 


0.18 


1.47 


11. 


45.00 


45.66 


2.35 


0.66 


12. 


46.00 


46.40 


1.31 


0.40 


13. 


47.00 


47.42 


0.96 


0.42 


14. 


49.00 


51.03 


0.12 


2.03 


15. 


50.50 


50.73 


0.50 


0.23 


16. 


51.00 


49.77 


0.31 


1.23 


17. 


51.00 


52.07 


1.29 


1.07 


18. 


51.00 


52.22 


0.69 


1.22 


19. 


53.00 


53.20 


0.72 


0.20 


20. 


57.00 


56.60 


1.22 


1.00 


21. 


63.00 _ 


63.34 


. 0.23 


0.34 


22. 


65.00 


62.20 


0.77 


2.80 


23. 


67.00 


64.98 


0.62 


2.02 



The exemplary data was computed for each breath contained in a 2-5 minute 
collection period. The calculated PaC0 2 values and standard deviations included in 
Table m were determined by analysis of the median six (6) values of ten (10) 
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consecutive Bfeaths. The standard deviations for the six (6) median values for each 
patient were computed for assurance of reproducibility of the data. The computed 
values were then compared to the PaC0 2 ABG values obtained simultaneously from an 
arterial blood gas sample. Using this method the accuracy for determination of 
5 PaC0 2 by the method of the instant invention is typically about ±2 mmHg. 

However, rather than outputting a single PeC0 2 value, the actual data range for a 
95 or 99 % confidence levels may be graphically displayed, as on the data display 
(Figure 3, item 11). 

Since other modifications and changes varied to fit particular operating 
10 requirements and environments will be apparent to those skilled in the art, the 
invention is not considered limited to the example chosen for the purposes of 
disclosure, and covers all changes and modifications which do not constitute 
departures from the true spirit and scope of this invention. 
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CLAIMS 



What is claimed is : 

1. A non-invasive arterial gas partial pressure determination system, 
comprising: 

expiratory gas sampling means for taking one or more samples of expiratory gases 
from a patient and developing a gas profile for at least one constituent of 
said sampled gas; and 

gas sampling data processing means for receiving data input representative of said 
gas profile, said data and selectively processing said inputed data to 
determine at least one arterial gas partial pressure value. 

2. The system of claim 1 , and further comprising means for assessing 
said gas partial pressure with respect to expiratory volume for said sampled gas. 

3. The system of claim 2, and further comprising means for evaluating 
said sampled gas in terms of gas partial pressure and expiratory volume iteratively 
during an exhale cycle to develop a multi-variable relationship. 



4. The system of claim 3, wherein said multi-variable relationship is 
20 developed in a time domain, further comprising means for converting said multi- 
variable relationship from the time domain to a volume domain. 

5. The system of claim 3, and further comprising memory means for 
storing a matrix of linear arguments recallable to form a vector, used by said data 

25 processing means in conjunction with said multi variable relationship to determine a 
partial pressure for arterial C0 2 . 



6. In a computer controlled data collection and processing method for 
monitoring arterial gas values for at least one gas, comprising the steps of: 
30 sampling a patient's exhaled breath to determine the partial pressure of said arterial 
gas in said exhaled breath as a function of incremental breath volume during 
an exhale cycle; 
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segregating said partial pressure data into three distinct phases representative of gas 

concentration of said at least one gas during said exhale cycle; 
extracting functional variables from at least two of said distinct phases in terms of 

slope and intercept values; 
recalling from memory a matrix of stored weighting parameters associated with said 

extracted variables; 
creating a vector based on said extracted variables modified by said stored 

weighting parameters; and 
determining an arterial gas value for said at least one gas based on said vector. 

7. The method of claim 6, and further comprising calculating said 
matrix of weighting parameters by an artificial intelligence regimen assigning weight 
to individual variables based on adduced significance. 



15 8. The method of claim 7, and further comprising sampling said exhaled 

breath at a frequency approximately between 10 and 1000Hz. 

9. The method of claim 8, and further comprising eliminating extreme 
calculated arterial gas values, statistically qualifying the average of multiple 

20 calculated preliminary arterial gas values, and determining the arterial gas value for 
said at least one gas. 

10. The method of claim 9, and sampling said exhaled breath at intervals 
corresponding to the requirements of the patient. 

25 

11. A system for deriving blood gas concentration for at least one gas, 
comprising: 

a) data acquisition means for measuring exhaled breath and 
determining concentration of one or more gas fractions thereof at volumetric 

30 increments of said exhaled breath; and 

b) data processing means for receiving said acquired data and 
extracting select parameters therefrom to form a relationship, and implementing said 
relationship to derive said blood gas concentration for said at least one gas. 
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12. The system of claim 11, and further comprising means for expressing 
the derived blood gas concentration as a partial pressure. 

13. The system of claim 12 wherein said at least one gas is C0 2 . 

5 

14. The system of claim 11 wherein said data processing means further 
includes a memory means for storing a matrix of weighting values to be applied to 
said relationship. 

10 15. The system of claim 14 wherein said weighting values are calculated 

by iterative assessment of predictive gas concentrations compared to actual gas 
concentrations for said at least one gas. 

16. A non-invasive method of determining blood gas content based on 
15 volume domain, comprising the steps of: 

a) taking time domain measurements of a plurality of a patient's 
expiratory breath during at least one respiratory cycle using a pneumotachometer, 
and a gas sensor for sensing the concentration of said blood gas in said breath, said 
measurements being at a frequency approximately between 10 and 1000 Hz; 
20 b) converting said time domain measurements to the volume 

domain; - 

c) mapping said volume domain data to form a curve; 

d) determining a threshold point for said concentration of said 
blood gas in said breath, whereat said curve transcends from a value below a 

25 threshold value to a value above said threshold value; 

e) calculating a first phase slope for said curve wherein said 
threshold point is an initial phase point, plotting subsequent points until said points 
deviate from a linear path, said deviation indicating the termination of said first 
phase; 

30 f) calculating a second phase slope for said curve by taking the 

last data point of expiration during said at least one respiratory cycle and regressing 
toward the termination point of said first phase; 
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g) extracting first phase and second phase slope values to create 
an expiratory relationship; 

h) combining said expiratory relationship with stored arguments 

into at least one vector; 
5 i) using said at least one vector to produce an estimate of the 

concentration of said blood gas in said blood; and 

j) visually displaying said gas content in said blood. 



10 



17. The method of claim 16 wherein said gas is C0 2 . 

18. The method of claim 17 wherein said threshold value is about 0.5% 
concentration of said blood gas. 

19. The method of claim 17, and further determining the partial pressure 

ofC0 2 . 

20. A non-invasive method of determining a dissolved gas content of the 
arterial portion of the pulmonary vasculature expressed as a partial pressure 
comprising the steps of: 

a) taking time domain measurements of a plurality of an adult 
patient's expiratory breaths using a pneumotachometer and gas sensor, at a 
frequency of about 100 Hz; 

b) converting said time domain measurements to a volume 

domain; 

c) mapping a curvature relationship of gas partial pressure versus 

gas volume; 

d) determining a threshold point, said threshold point being the 
point where said cuivature^lationship : transcends from a value below 0.5% of said 
gas partial pressure to a value above 0.5 % ; 

e) calculating a first phase slope wherein said threshold point is 
an initial point including plotting subsequent points until said points deviate from a 
linear path, said deviation indicating the termination of said first phase; 
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f) calculating a second phase slope by taking the last data point 
of expiration and regressing toward the termination point of said first phase; 

g) using the first and second phase slopes to quantify a expiratory 
gas relationship; 

5 h) combining said relationship with stored arguments into at least 

one vector; 

i) using said at least one vector to estimate the partial pressure 
of said arterial gas; and 

j) visually displaying said partial pressure of said dissolved 

10 arterial gas. 



21. 



The method of claim 20 wherein said dissolved arterial gas is C0 2 . 
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